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Abstract 

A new mean-field type theory is proposed to study order-disorder transitions 
(ODT) in block copolymers. The theory applies to both the weak segregation 
(WSL) and the strong segregation (SSL) regimes. A new energy functional is 
proposed without appealing to the random phase approximation (RPA). We find 
new terms unaccounted for within RPA. We work out in detail transitions to the 
lamellar state and compare the method to other existing theories of ODT and 
numerical simulations. We find good agreements with recent experimental results 
and predict that the intermediate segregation regime may have more than one 
scaling behavior. 



PACS numbers: 83.70.Hq, 36.20.-r, 61.25.Hq 



* depablo@engr.wisc.edu, rebei@cae.wisc.edu. 



1 



1 Introduction 



Predicting morphologies of block copolymers continues to be a very challenging problem 
from the computational point of view. To make this latter task more feasible 'nice' energy 
functionals are needed and this is the goal of this paper. For linear copolymers, the self- 
consistent field method (SCFT) developed by Helfand and others is considered to 
be the method of choice. It applies quite well to all regimes of segregations. Recently it 
was successfully used to predict different morphologies of linear triblock copolymers ||. 
The method is based on a naive mean field approximation to the partition function. 
Fluctuation effects can also be added but at the expense of complicating the method ||. 

Around the order- disorder transition point, Leibler's || field theory for diblock 
copolymers is superior to the self-consistent field theory. For this reason there have 
been attempts to generalize it to all types of segregations, intermediate as well as strong 
segregations. Ohta and Kawasaki |T0|-p^1 were the first to propose a free energy func- 
tional of Leibler's form that treats the strong segregation case of a diblock copolymer. It 
gives comparable results as the self-consistent method and its predictions are well sup- 



ported by experiment [pL3|, pL4|| - Fluctuation effects were later added to Leibler's theory 
by Fredrickson, Helfand and Barrat [|I5|--|T9|| who showed that the peak of the scattering 
function not only depends on %iV, the Flory-Huggins parameter, but also on the average 
segment length and volume. The Leibler free energy functional was also shown to be 
useful in the strong segregation case if the wave vector dependence of the energy func- 
tional is fully kept , since higher order spatial harmonics of the order parameter become 
more and more important as the temperature is lowered in the ordered phase [pO| .|2T|. 

Here, we set to find a similar, but simpler, energy functional for incompressible 
diblock copolymers that is valid for both the weak and strong segregation regimes. Our 
work has the same spirit as the recent work of Stepanow |22| where he used graphical 
methods to obtain an improved self-consistent expression for the structure factor of a 
diblock copolymer. He avoided using the random phase approximation (RPA) [^3| , i.e., 
Leibler's theory, and instead finds an expansion of the partition function in terms of an 
effective potential p4[-f26|1 . Here we adopt the same goal of developing an expression of 
the free energy based also on an effective potential. The random phase approximation 
will not be used to get our energy functional. 

Besides succeeding in finding a new expression for the energy of a diblock copolymer, 
we also find new "ideal " terms missed by RPA and already found by Holyst and Vilgis 



27fl in the polymer mixture case. 

Our method is strictly functional; we do not use any graphs. Our results to lowest 
order are similar to Stepanow result [2l| and our energy functional is qualitatively the 
same as Leibler's. The quartic term is similar to Leibler's fourth order term but is much 
simpler to work with. 

The paper is organized as follows. In section II, we develop the formalism. The 
details are left to the Appendix. In section III, we solve the transition to the lamellar 
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morphology of a symmetric diblock. In section IV, we compare our results with those 
of the self-consistent method, Leibler's method, and simulations [p8 . 



2 Free Energy of a Block Copolymer 

We start by deriving a new expression for the free energy of a block copolymer melt. 
Even though in this paper we are mainly interested in incompressible diblock copolymers, 
it can be easily generalized to, e.g, triblocks. The method we use bypasses the random 
phase approximation and the use of virtual sources, as was done originally by Leibler |§. 



The Hamiltonian we use for our system is that of Edwards [24] generalized to copolymers 
and can be taken of the form: 



H = H + V, 

n pN 



*oN = ^£/ o dr (-11 j 

i r 2 

V[r,] = - / drdv' P,(r)^(r')V 2J (r,r') 



i,j = A,B (1) 

Ti is a curve along the i'th macromolecule that has two types of monomers A and B with 
densities Pa( y ) and Pb( y ), respectively. We assume both monomers to have the same 
Kuhn length o. There are n chains in the melt, each with f (f —fN) monomers of type 
A and (N — f) of type B. The interaction potential V is taken to have the simple form 



V(r,r') = Po ( ° x lyinis)-^')), 



(2) 



where \ is the Flory-Huggins constant and po is the total average density of monomers. 
The partition function of this incompressible system of macromolecules is then given by 



Z = j d(r t )6(l-p A (r)-p B (r))exp -\H + ^J drdr' 7? T (r)V(r, r')pV)]) 



(3) 



In the above we have set the Boltzmann constant to one and used a vector notation for 
the densities, i.e. 
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The densities are given by 



N n f f 



(5) 



and 



n n r 1 

^(0 = -E / ds5(r-r a (s)). (6) 

In Z, the second term in the exponential is a quadratic symmetric form and hence it 
can be diagonalized. This diagonalization allows us to deal with a new virtual set of 
monomers that are 'decoupled'. We therefore introduce a new set of variables p\ and p2 
such that: 

Pi( r ) = \{Pa{t) +p B (r)), (7) 

P2{r) = \{ PA {r)-p B {r)). (8) 
In this new set of variables the potential becomes diagonal with 

Vii = 2poX, (9) 

V22 = -2xPo, (10) 

V 12 = V21 = 0. (11) 

Ignoring the incompressibility condition for now and after isolating the free energy of 
the disordered state we have 



where 



and 



x exp{-±p° V^Apg}, (12) 



Vri = d(ri) exp(-ifo), 



A Pa (r) = p a (r) - p° a , a = 1,2. 
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For a symmetric diblock, 



(13) 



Upon introducing a two-component Hartree-type field tp a with which these new virtual 
molecules are interacting, the partition function becomes 



+ifdr[<p a (r)+ip°V af3 ]-Ap a (r)}. 



(14) 



We introduce now two more fields, $ and p. Since 

J V<$> a <S($ Q (r)-Ap Q (r)) = l, 

then we have 



(15) 



Z = exp(-F ) J VipV^VfiV( ri )exp{-ifi a ^ a } 
x exp-J-^V"^ + i(p a + p a + ip°pV al3 )Ap a }, 



(16) 



where 



F =n X Nf(l-f). (17) 
is the energy of the disordered state. Now, we make a change of variables and let 



(Pair) = ¥a{r) + p a (r)+i J dr l/ Q/3 (r)pJ(r). 



We also define a new functional 



ln^ D(ri)exp -^H (ri)+i J dr ip(r) ■ p(r) 



(18) 



(19) 



This functional can be expanded in around a homogeneous state 



F M = E™=o ^ E M //•••/ dr 1 dr 2 ...dv m C ala2 ... am (r 1: r 2 , .., r m ) 

0.2 (^2) ■■■(Pa.m (j^m) 
OLi = 1,2. 



(20) 
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The coefficients C aj g„. are given in Appendix A. In the following we keep only terms up 
to the fourth order. After some rearrangements and writing ip in place of Tp, we have 
the following expression for the partition function, 

Z = J V^VnV<pexp{~<p a {V^ + C a p)ipp (21) 

1 r, - 1 

x J exp{^C a/3A7 (^ a;3 - ip a ipp)(ty X ~/ ~ f\f-y)}- 

Here we have introduced a new pairing field ty a p so we can cancel the quartic term in 
<p a . We have also set 

iE(v) = ip° a + Vjp (r), (22) 

and have not written the space integrals explicitly. Then, we can integrate the (p field 
and then the p field exactly. We are left with only two fields $ and \& ( see Appendix 
B) 

Z = J V$Vtyexp{-F{<$>,^)}, (23) 



where 



+ \(«B-M + \ lo S det A <*P + ^ lo § det ( 24 ) 



and 



C«(r) = S (r) + p° - IC/^aVa, (25) 
A={l + \^U)U-\ 



b = v^ 1 - r^'V, 

U a /3 = (V a [3 + C a p) . 

In the above A, B, A, and U are spatial dependent matrices. U a p is the desired effective 
potential ||24j| . We expand the logarithmic terms in powers of U. We do the same when 



we seek an expression for B. In all expansions, we keep only quadratic terms in . The 
final lowest order expression we get from this expansion is our expression for the energy 
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"<•) = 2^E*w(z4)- 2wo) * ( - ,) (26) 



1 y^^/ ,C 2 222(q,-q,p,-p) 



3!M) 2 ^ w [C 22 (g)] 2 C 22 (p) 
1 1 y-v C' 2222 (gi, g 2 , g 3i -qi-q 2 - qs) 



4! (W) 3 ^ C , 22 (gi)C 22 (g 2 )C 22 (g 3 )C 22 (gi + g 2 + g 3 ) 
x$(gi)$(g 2 )$(g 3 )$(gi + g 2 + g 3 ). 

It can be easily shown that wo//(2nC 22 (g)) is actually Leibler's structure function. 
Leibler's structure function S(q) is given by (Appendix A) 

n-i ( s = S AA (q) + 2S AB (q) + S BB (q) 

{q> S AA (q)S BB (q) - S AB (q) ' 1 ' 

This is indeed equivalent to C^f(q) with the proper normalization. However it is easier 
to see this graphically and in fig.l, we plot this function. Hence the first and the second 
terms are the usual RPA result for the inverse scattering function. The third term is 
new and is not captured by the RPA approximation. It is not due to fluctuations. It 



has been first pointed out by Holyst and Vilgis in their study of polymer blends |27 
The fourth term is familiar but the coefficient is much simpler to calculate and behaves 
differently than Leibler's term for large wave vectors. In fig. 2, we plot both terms for a 
subset of wavenumbers. Our fourth order coefficient is given by 
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Figure 2: Comparison of fourth order coefficients of this theory and Leibler's RPA 
theory. 



9 



nU)/ \ C 2 222(ll, 12, g 3 , ~ gl — g2 ~ Q3) / oc n 

C ( ; (gi, g 2 ,g3, - gi - ga - ga) = ^ ( \n ( \n ( \n ( — 7 i V ( 28 ) 

C22(gl)C 2 2(g2)C 2 2(g3)C 7 22(gl + <?2 + g 3 ) 

This term is independent of any three-body correlation functions. Leibler's I^ 4 ) term 
does however depend on three body correlation terms, and this greatly complicates 
computations involving it; it is given by |9j 

T (4) (gi, g2, ga, - gi - g2 - ga) = lijki(qi, ga, g3, - gi - g2 - g3)(^ x (gi) 

-^(gi))(^i(g 2 ) - S^fc)) x (^(g 3 ) (29) 

- 5 feife))(5 M 1 (g4) - S^fa)). 

with k,l = A, B. ^%ju is a function of two-point, three-point and four-point correla- 
tion functions. In fig. 2, we plot both coefficients for gi = — g 2 = g3 = — g4- 

To lowest order, the propagator of this theory also agrees with that of Stepanow |2"2" 



except that our self-energy term has additional contributions from the fourth term, 
Tr(AU), of the free energy in Eq.(B-7). Hence to consider any fluctuations this term 
must be included from the outset and therefore the work in [2^] is more thorough than 



the one presented, e.g., in |[J| and others based on RPA. In the next section where we 
use our functional to study transitions from a disordered state to a lamellar state, we 
will not include the quadratic non-RPA term in order to compare our functional with 
that of Leibler's. 



3 The Lamellar Solution 

Following Melenkevitz and Muthukumar [EDI, we conjecture a solution that minimizes 



the energy functional. Knowing that in SCFT the densities are found by solving a 
modified heat equation, we choose a function of the following form 

^ 2 1 

®( x )= Yl ^ ex p[-2^ A ) 2 ] sin ^ x )- ( 30 ) 

1=1,3,... 

This choice is also dictated by the fact that 

dx${x) = 0, (31) 



and the solution must be periodic. The wave vector q is given by 

e = ™- < 32 > 

where D is the lamellar periodicity. A is another parameter besides D that is related 
to the wall thickness of the interface region between the two components of the diblock 
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copolymer. Both parameters are to be found by minimizing the energy H with respect 
to them. Actually, we solve for A for a given D and then calculate the corresponding 
energy and choose the solution with the lowest energy. We rescale dimensions in terms 
of the radius of gyration so the energy per chain is given by 

b m b p b r b s C 2 222(m,p, r, s) 



384 m »fri 3 ... c Am)C^(p)C^(r)C^(a) 



P 

r 

-m—p—r 



where 

b m = A exp [_i( 27rm A/ J D) 2 ]. (34) 

Figures 3-8 summarize all the results about this particular solution. In particular, 
we observe that this solution predicts that the order- disorder transition (ODT) occurs 
for xiV right below 10.5. Immediately below the transition temperature, our energy 
functional shows that the behavior of the polymer chains is no longer Gaussian. For 
Gaussian chains, the scaling factor 5 between hi(D) and ln(xiV) is zero. In our case, we 
find that for xiV less than 13, 5 is approximately equal to 0.26. For xiV between 13 and 
30, 5 = 0.52, and for \N above 30, 5 becomes about 0.19. Hence the 2/3-power law 
between D and N in the SSL is also verified by this solution. So far only one intermediate 



region has been observed ||13|| . Here, this solution suggests that the intermediate region is 
really more than one. Since the behavior below the ODT is believed to be nonuniversal, 
it will be interesting to see if this predicted behavior is also observed in some symmetric 



diblock copolymers other than the one treated in JT^]. The inclusion of the non-RPA 



term that we omitted in this solution will not change this overall picture. It has only a 
moderate effect at high values of \N where higher and higher wavenumbers are needed 
in the energy. 
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Figure 3: Density profiles for the lamellar morphology given by Eq.(30) for \N = 
10.5,11,12,25,50,100. 
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4 Comparison and discussion: 



In fig. 3, we plot density profiles for different X-^' s - We observe that the order- disorder 
transition occurs below \N = 10.5. The segregation amplitude grows much faster for 
xN just above the ODT temperature x^ c and less than 30. Hence in this theory, 
the strong segregation regime is attained much faster than the self-consistent method 
(SCFT) predicts This is also confirmed by recent experiements on symmetric diblock 
copolymers [O]. Figure 4 compares our results to that of the SCFT calculations at 



xN = 12.5 and simulations Pq| . The simulation was done with chains of forty-eight 
segments. Assuming that the critical temperature, xN c , for these chains is also close to 
10.495 as is the case for infinite chains, we find that the simulation gives a result that 
falls between our result and the SCFT result. However the ODT temperature for these 
short chains is expected to be less than in the ideal case. The SCFT curves were found 
by solving Eqs.(3-7), that appear in 0, using a finite difference method and using a 
random configuration as an input. Our results agree well with those given in 0. In 
fig. 5, we instead predict the xN value for which our theory coincides with the simulation 
results. We find a xN value of approximately 11.15. Assuming now that x^/x^c — 1-2 
corresponds to xN = 11.15 in our theory, we plot in fig. 6 and fig. 7 density profiles for 
xN = 22.3 and xN = 44.5 and compare our results against those of simulation. We 
find relatively good agreements between theory and simulation especially far from the 
interfaces. Finally, in fig.8, we check if our theory predicts the observed scaling behavior. 
Clearly, we can distinguish three regimes from the plot. The intermediate regime extends 
from about xN = 13 and to about xN = 27, has a scaling factor 5 of about 0.52. 



Above xN = 30, the scaling factor is about 0.19, in agreement with observations |T3 
Compared to Leibler's energy functional [^TJ, our results are much closer to the SCF. 
This is due to the fourth order term which is smaller than Leibler's fourth order term 
for large wavenumbers. For large xiV, our energy functional gives closer periods to the 



SCFT method than the full Leibler Hamiltonian [21 1 that overestimates periods by as 



much as 30 %. Moreover, our theory predicts that transition to the strong segregation 



regime occurs at xN ~ 30. In complete agreement with the experimental results in [13 



In the SCFT, the SSL is believed to occur around xN = 50. Ref. |2tJ suggests that the 
SSL starts to occur for %iV larger than 90. 
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Figure 4: Density profiles of the symmetric lamellar solution for xiV = 12.5. The 
Simulation is done for xiV/xiV c = 1.2. 
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Figure 8: lot of InD vs. ln(xiV) for a symmetric lamellar morphology. For high xN, the 
scaling factor is approximately 0.69. 5 is the slope of the best fit segment to the data in 
three different regions. The weak segregation limit (WSL), the intermediate segregation 
configuration (ISC) and the strong segregation limit (SSL). 
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Appendix 
Appendix A 

Let Q be a partition function of a melt of polymer chains in external fields <pi and <pi 

Q[(pi] = J V(r) exp{-^iS0 - iy? 2 A0} (A-l) 

where 

A</> = (1 - - f<p 2 (A-2) 

</?i, v?2 are densities. 

We next expand logQ in terms of <fi and <p 2 ■ Since Q should be invariant under 
A0 — > —A(f) and £0 — > — £0, only even powers of </?j are present, hence we write: 



log<5 = ^ / dxidx 2 C a p(x 1 ,x 2 )(p a (x)(pp(y) 



(A-3) 



+ 



+ 



-jj y dXi...dx 4 C a( g 7 5(xi, ...,l4)^(ll)...^(l4) 



where, e.g., 



C 22 (x 1 ,x 2 ) =< A(p(x 1 )A(p(x 2 ) > 

C 2 222(^l--^4) =< Ay>(xi)...Ay>(x 4 ) > 

- C , 22 (a:i,X2)C , 22(a;3,^4) 

- C 22 (£i,Z 3 )C220e2,£4) 

- ^22(^1, X4)C 2 2(X 2 ,X 4 ) 



(A-4) 
(A-5) 



The averages < ... >o are evaluated with a Gaussian distribution. Hence we find 

C 22 (q) = S AA (q) - 2S AB (q) + S BB (q) 

with 



Saa(x) = -^(f A (x) + exp(-f A (x) - 1) 



S AB (x) 



x 2 
-1.0 



(A-6) 



x z 



(exp(-f B (x)) - 1 - exp(-x) + exp(-f A (x)) 



and x = q' 



,a 2 N 
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Similar expressions for the other correlation functions, such as C2222? can straightfor- 
wardly derived, but with much more labor. It amounts to calculating all correlation 
functions, G a/ 3 1 s(xi,X2,x 3 ,X4 : ), in the Gaussian distribution of the form 



Gap-ys(xi,X2,x 3 ,X4) = J ds 1;a J ds 2 ,/3 J ds 3j7 J ds4,s{S(xi - x(s 1>a ))S(x 2 - (A-7) 
x(s 2 ,/3))S(x 3 - x(s 3j7 ))5(a;4 - x(s 4 ,s)))o 
where a, (3, 7, 5 — 1, 2. Hence C2222 will be a linear combination of all these functions. 



Appendix B 

We start from 



= J V$Vm exp{-ip a <S> a - ^p a Vjp p - ip° a p a } (B-l) 
x J Dv ? exp{-iv9 a (V^ 1 + C a p)ifp + iZ a (p a 



The integral over if is Gaussian and can be done exactly. If we set 

1 
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A a /3 = V a p + C a/3 + -^xsC a/3 x5 (B-2) 



we have: 



Z= j V<S>VV exp{^ a pC a p X ^ x , - \ P lA-y p - |Tr log - A aP } 

x j Vp exp{-\fi a V~^np - ip a $ a - ip° a p a + lPaV~ (3 1 A^V x ~ 1 p 1 (B-3) 

+^v-^a- p IpI} 

Again the p— integral is Gaussian. After integrating out p, we have 

Z = J P$P*exp{^* a/3 C^ A7 * A7 (B-4) 
l /*A- a lpl- l -Tr\ogA aP 



If we set 



T- 1 = V- 1 - V^UV- 1 , (B-5) 
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leave out terms of order (C Qi g>,,) 2 , and use the incompressibility constraint, i.e., 

P°2 = (B-6) 
$i = 

we get after some quite heavy algebra the following simple form: 

.F($, A) = l -J d Xl dx 2 <S>{l)T{l, 2)$(2) (B-7) 

dx 1 ...^ 4 A(l,2)C 22 1 22 (l,2,3,4)A(3,4) 

- X - J dx 1 ...dx 8 ^(l)(TV- 1 U)(l, 3)A(4, 5)(f/\/~ 1 T)(5, 8)$(8) 

+ ^/ dx 1 dx 2 A(l,2)U(2,l) 

l -j dx 1 ...dx 6 (T\/- 1 f/)(l,4)A(4,5)(f/\/- 1 T)(6,l) 



+ 12 



where 



and 



A(xi,x 2 ) = J dx 3 dx 4 C 2222 (l, 2, 3,4)^(3,4) (B-8) 



= $ 2 (:r) 

Now we can integrate over A since the integral is only Gaussian, and we find that 

Z = J V$exp{-F($)} (B-9) 



where 



A.F($) = - F (B-10) 



2voZ 



1 1 >p C 2222 (g, -p) v 
6t;o/ Z. [C 22 (g)pc 22 (p) 1 [ q) 

l l v-^ C 2222 (q,p, fc, -(<? + p + A;)) 



E 



and 



24 H/) 3 ^ C 22 (g)C 22 (p)C 22 (A;)C 22 (g + p + k) 
x $(g)$(p)$(/c)$(-p -q-k) 



T{q) = c4) + M,) (B ' U) 



is now the effective potential for the two-component incompressible copolymer melt. 
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